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Abstract 

We study the effects of tlie MSSM contribution on i? — > (j)K decays using the perturbative 
QCD approach. In this approach, strong phases can be calculated, so that we can predict the 
values of CP asymmetries with the MSSM contribution. We predict a large relative strong phase 
between the penguin amplitude and the chromomagnetic penguin amplitude. If there is a new CP 
violating phase in the chromomagnetic penguin amplitude, then the CP asymmetries may change 
significantly from the SM prediction. We parametrize the new physics contributions that appear 
in the Wilson coefficients. We maximize the new physics parameters up to the point where it is 
limited by experimental constraints. In the case of the LR insertion, we find that the direct CP 
asymmetries can reach about 85% and the indirect CP asymmetry can reach about —30%. 
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I. INTRODUCTION 
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B — s> (j)Ks decay may be useful in the search for new physics beyond the standard model 
(SM). The time-dependent CP asymmetry of decay into CP eigenstates fcp can be 
written as aj^pit) = Aj^p cos^AMst) + Sj^p sin(AMBt), where Af^p and Sj^p characterize 
direct CP violation and indirect CP violation, respectively. In the SM, both Aj/^Ks 
A(f,Ks vanish, and both Sj/^Ks S^Ks must equal to sin(20i). Any difference between 
Sj/ipKs ^<l>Ks larger than 0(1%) would be a signal for physics beyond the SM l|. The 
B (pKs decay amplitude is induced only at the one-loop level, so that new physics might 
contribute to this decay through quantum effects. At present, BaBar and Belle collaborations 
have the following results 0, [s | : 

' 0.741 ±0.067 ±0.034 (BaBar, 81 fb"^), 
0.733 ±0.057 ±0.028 (Belle, 140 fb"^) , 
with Aj/^Ks = 0- lu the (pKs mode, they have reported the following results 

0.38 ±0.37 ±0.12 (BaBar, 110 fb"^) , 
-0.15 ±0.29 ±0.07 (Belle, 140 fb"^) . 

0.45 ±0.43 ±0.07 (BaBar, 110 fb"^) , 
-0.96 ± 0.50 IH? (Belle, 140fb~^) . 
Belle collaboration found a 3.5(T deviation from the SM prediction. In contrast with Belle, 
the BaBar result is consistent with the SM. If the Belle result continues to hold, and both 
experiments agree, then their result is a signal for new physics. There have been many 
papers which studied new physics contributions to i? ^ (pKs decay amplitude. Some 
authors analyzed the supersymmetric contribution using the mass insertion approx- 
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imation, which is a powerful tool for model-independent ana. 
with the minimal supersymmetric standard model (MSSM) 



/■sis of new physics associated 
j]. The new physics contribu- 



tions come into the Wilson coefficients, which can be calculated perturbatively The 
problem is how to calculate the decay amplitudes with nonperturbative contributions. To 
calculate the decay amplitudes, some authors used naive factorization jll|, generalized fac- 



torization 



12i |. or QCD factorization [13]. Each method is plagued with large theoretical 



uncertainties. There are other approaches for this calculation, one of them is the perturba- 
tive QCD (PQCD) approach [1^. In this paper, we use the PQCD approach and the mass 
insertion approximation for an estimation of the MSSM contribution in i? ^ (f)K decays. 
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FIG. 1: These diagrams generate a large strong phase in the PQCD approach. The dashed hnes 
denote the physical intermediate states. 

The PQCD approach for exclusive B meson decays is based on the factorization theo- 
rem, in which the decay amphtudes can be separated into perturbative and non-perturbative 
parts The non-perturbative parts are factorized into meson wave functions, which are 
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derived from the other methods, for example, the light-cone QCD sum rules 
strong phase, which comes from physical intermediate states, is calculable in the PQCD ap- 
proach. A large strong phase is induced from an annihilation diagram such as Fig. H^a) . 
The strong phase mainly comes from the cut on the virtual gluon line. The source of the 
strong phase is one of the important differences between the PQCD approach and other 
methods. PQCD has been applied to some hadronic two-body B meson decays at lead- 
ing order in a,, and the results are consistent with experimental data except for rj^'^K 
and 7r°7r° \lA |l9|. B — > (j)K decays were also calculated using the PQCD approach in 
Refs. QH- 

Another important difference between PQCD and other methods is how to calculate 
magnetic penguin diagrams that are induced from the chromomagnetic penguin operator 
OgG. In many models, the chromomagnetic penguin amplitude is most sensitive to new 
physics 0, 23 1 . However, it is difficult to calculate the chromomagnetic penguin in naive 
factorization and generalized factorization, because we do not know the magnitude of g^, 
which is the momentum transferred by the gluon in the chromomagnetic penguin operator. 
Therefore, proponents of these factorization approach treat as an input parameter, so that 
the result is directly proportional to the assumed values for (l/q^) 'l2|. In the PQCD and 
QCD factorization approaches, the chromomagnetic penguin amplitudes can be calculated 
without any assumption for the value of q^. We studied the chromomagnetic penguin using 



the PQCD approach in Ref. 



24| , and we found that the chromomagnetic penguin generated 



a strong phase from the diagram as Fig. db). In the PQCD approach, is written as 
(1 — X2)x^M'^ — \k2T — fcsTp- Here, X2 and X3 are momentum fractions of partons in K 
and mesons, respectively. fe2T and k-^T are transverse momenta of the partons. In QCD 
factorization, can be written in terms of the momentum fraction of partons too, however, 
they expand the amplitude in power of \k2T — k-iT\'^/[{^ — X2)x^M'^ and for the leading order 
they have = (1 — X2)x^Mq. In this expansion, never vanishes. There is no absorptive 
part in the amplitude and the strong phase is not generated from the chromomagnetic 
penguin. The fact that we get large imaginary part implies that this expansion is not valid. 

The outline of this paper is as follows. First, we consider the MSSM contribution in 
the effective Hamiltonian for B meson decays. We present the Wilson coefficients with the 
MSSM using the mass insertion approximation. Next, we briefly review the PQCD approach 
for the exclusive B meson decays. We show the result at the leading order in Us and how to 
calculate the chromomagnetic penguin amplitudes. Furthermore, we calculate the MSSM 
contribution 'm. B (f)K decays with the LR, RL, LL, and RR insertions. We consider the 
both B^ (pKs and B^ —>■ (pK^ modes, and calculate the branching ratios, the direct CP 
asymmetries, and the indirect CP asymmetry. Finally, we summarize this study. 



II. MSSM CONTRIBUTIONS m B ^ (f>K DECAYS 



A. Effective Hamiltonian for B Meson Decays 



We use the effective Hamiltonian in the calculation of B meson decays j25[. The Hamil- 
tonian is expressed as the convolution of local operators and the Wilson coefficients. The 
effective Hamiltonian with a 6 — s transition is given by 
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E V*,,V,.,(CM0^^'\^^) + C2{^^W 
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\i=3 



+ H.C. , 



(4) 



where V*,^ and Vq'b are the Cabibbo-Kobayashi-Maskawa matrix elements [20]. Oi_io are 
local four-fermi operators, Oj^ is the photomagnetic penguin operator, and Osg is the chro- 
momagnetic penguin operator. The local operators are given by 



or 



(9') 



{siq'Av-A{q'ibi)v- 



{siq'i}v-A{q'jbj)v 



Oz = {Sihi)v-A Y,iQj<lj)v-A , O4, = {Sibj)v-A Y,i<lj<li)v-A , 

q q 

O5 = {sibi)v-A^{qjqj)v+A , Oq = {sibj)v-Aj2i^j^i)v+^ ' 
q q 
3 3 
O7 = -isibi)v-Aj2eq{qjqj)v+A , Og = -{Sibj)v-Aj2^q(^j'ii)y+A ' 
q q 

3 3 
O9 = -{Sibi)v-Aj2^q(^j'ij)v-A , Oio = -{sibj)v-A^eq{qjqi)v-A , 

^ q ^9 

O77 = 8^^6^i^'''(l + l5)biF^u , OsG = -^m5S,(7'''^(l + l^)T^bpi^ . (5) 

Here, i and j are color indices, q is taken to be d, s, and c, and (gg)y±A = ?7/i(l i Ts)?- 
We define the covariant derivative as D^j, — + igsT"-A^ — ieA^j^, so that the signs of 
the magnetic penguin operators are different between O-j^ and O^g- We integrate out the 
degree of freedom of high energy particles, then the Wilson coefficients include high energy 
information. If we consider the new physics effect on B decays, then we need to calculate 
the Wilson coefficients with new physics contributions. 



B. MSSM Contribution and Mass Insertion Approximation 

We consider the MSSM contribution for B — > (t)K decays. Generally, there are new 
sources of the CP violation and the flavor changing neutral current (FCNC) in the MSSM, 
so that there may be direct CP violations, and it is possible that S^Ks becomes different from 
Sj/^Ks- Since we do not want our computation to depend on specific SUSY models, we use 
the mass insertion approximation to calculate the Wilson coefficients with the MSSM. In the 
mass insertion approximation, the FCNC effect appears in the squark propagators through 
the off-diagonal elements in the squark mass matrices. The decay amplitudes are expanded 
in terms of {S'^)ij — {Vjm^Vci)ij/m^, where m| is the squared down-type squark mass matrix, 
rriq is an average squark mass, and Va is the matrix which diagonalizes the down-type quark 
mass matrix. Of course, we must consider the region of < 1. For example, a transition 

of a right-handed fermion to a left-handed fermion is parameterized by There are 

four mass insertions: {S'lj^)ij, {Sjiji)ij, {^tidij-i {Sjii^)ij. The 6 — > s transition is induced 
from the gluino-squark loop, the chargino-squark loop, and the neutralino-squark loop. In 
this study, we consider only the gluino contribution, which is dominant in many models. 
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The Wilson coefficients for penguin and magnetic penguin are given by 
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at the first order in the mass insertion approximation 



(6) 



Here, Ms is the SUSY scale, and 



X = m'~/m'~, where rrig is the gluino mass. B{x), P{x) and M(x) axe the loop functions, which 
are calculated from box diagrams and penguin diagrams Q,|23|. New physics contributions 
induce additional operators, which are obtained from Eq. (jS)) by exchanging L and R. The 
penguin coefficients C^^{Ms) — Cq^{Ms) depend only on 5^^, and the magnetic-penguin 
coefficients Cf^{Ms) and C^qIMs) depend on both 6'[^ and 5^^. It is noted that the 
terms have a chiral enhancement factor mg/mb. This factor is important when we constraint 
the parameters from the branching ratio for B Xg'j, as we will return it on later. 



III. B ^(t)K DECAYS IN PQCD APPROACH 

Let us briefly review the PQCD approach for exclusive B meson decays. The PQCD 
formalism is based on the factorization of decay amplitudes into a product of long-distant 



physics, which is identified with meson wave functions, and short-distant physics [28[. The 
meson distribution amplitudes are universal in the processes under consideration, and they 
are determined from experiments and/or other theoretical methods: the light-cone QCD sum 
rules, lattice calculations, etc. Process dependence is exhibited in the short-distant part. 



^ The signs of {Sf 1^)23 terms are opposite from those in Ref. Pi and the same as those in Ref. [lo| . 
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FIG. 2: A leading-order diagram in B ^ (j)K. 

For example, let us consider B (pK. Figure. |21 shows the PQCD quark-level diagram for 
this decay. At the black blob, the decay b sss takes place. The ss pair forms a (p meson. 
The other s quark recoils against the meson carrying almost m;,/2 momentum. In order 
for the spectator quark to form a K meson together with the fast moving s quark, it has to 
exchange a gluon in order to change its momentum from ki to /c2- 

We consider the B (pK decays. In the light-cone coordinates, the B meson momentum 
Pi, the K meson momentum P2, and the meson momentum P3 are taken to be 

Mb, , Mb, o , Mb , o 

Pi = -|(1,1,0t), P2 = ^(l-r5,0,OT), P3 = ^(r5,l,OT), (7) 

where = M^/Mb- Here, we consider the B meson to be at rest, and the K meson mass 
is ignored. The momenta of partons ki, k2, and ks defined in Fig. |21are written as 

ki = (0, xiPf , kir) , k2 = (X2P2+, 0, k2T) , ks = (0, xgPg", ksr) • (8) 

A hard part in Fig. |21 has two propagators: 

1 1 1 1 , , 

{ki-k2f (Pi-A;2)2-M| -X1X2MI- |fciT-fc2TP -X2MI - \k2T\'' ' ^' 

If we neglect transverse momenta of partons, then a singularity arises from the end-point 
region of parton momenta since the hard part is proportional to l/xix\. Therefore, we con- 
clude that the transverse components are present and they cannot be ignored. Retaining kr 
for partons, large double logarithms appear through radiative corrections. The resummation 
of those double logarithms leads to the Sudakov factor 2^. The Sudakov factor suppresses 
the end-point of the parton momentum and the large transverse separation between a quark 
and an antiquark in mesons. Therefore, it guarantees a perturbative calculation of the hard 

r~i 

part [30]. The other double logarithms appear from the end-point region of the parton mo- 
menta. Their resummation, which is the so-called threshold resummation, leads to another 
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FIG. 3: Leading order diagrams for B — > (pK decays. The ellipses denote the meson wave functions, 
and the thick points denote the local operators. 



factor in the hard part [31]. This factor also ensures the absence of the end-point singular- 
ities in PQCD, and the arbitrary cutoffs used in QCD factorization are not necessary. A 
typical decay amplitude for B —>■ (pK can be expressed as the convolutions of a hard part, 
meson wave functions, and the Wilson coefficient, in the space of Xj and 6j, where 6, is the 
conjugate variable to kiT j3- 



M = [ dxidx2dx3 [ dbidb2db3^K{x2,b2) e •^^^^^'^^'^^ $<i(x3, 63) e 
JO Jo 



xC (t) H{xi, X2, X3, 61, 62, bs, t) J(xi, X2, X3) 61) e-^^("i'''i'*) . (10) 

Here, and $0 are the meson wave functions, and H is the hard part. Sb, Sk, 

and denote the Sudakov factors, and J denotes the threshold factor. The scale t, which 

Mb — rrify. The formula 



characterizes the hard part, is of order of y AMb where A 
in Eq. (jlOj) is the typical amplitude for an exclusive B meson decays based on the fcr 
factorization. 



A. B ^ (j)K in the Standard Model 



For B —>■ (j)K decays, we calculate the diagrams shown in Fig. 01 at leading order in the 
PQCD approach. The diagrams (a) and (b) are dominant contributions, and the diagrams 
(e) and (f) generate a large strong phase 1^. Except for a small tree diagram contribution 
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in the (pK^ decay amplitude, both 5° and (pK^ decay amphtudes get 

contributions from pure penguin graphs. For B meson, we use the model wave function 



(l)B{x,b) = Nbx'^{1 — xY exp 





(xMb^ 




H 




1 2 



(11) 



where h is the conjugate space of fcr- ^b is the shape parameter to be 0.40 ± 0.04 and A''^ 
is the normalization constant 1331. For K and 6 mesons, we use the wave functions that 

amplitudes for B — > (pK are shown in Ref. [20|^. We get the following numerical results 
within the SM: 

Br(S° ^ ^K^) = (8.5l|[! ± 2.6) x 10"^ , (12) 
Bt{B^ (f)K^) = (9.3l|} ± 2.8) X 10~^ . (13) 

The values for various parameters that we use in this calculation are presented in Ap- 
pendix El Here, the first error is estimated from the shape parameter ub in the B meson 
wave function, and the second error comes from higher-order contributions. We expect that 
the higher-order contributions are about 30%. The theoretical errors are reduced in CP 
asymmetries, because errors associated with wave functions cancel out between the denom- 
inator and the numerator. The current experimental data are given by 

I (8.4+H ± 0.5) X 10"^ (BaBar) , 
y (9.0l?J ± 0.7) X 10-6 (Belle) , 

f (10.0+H ± 0.5) X 10-6 (BaBar) , ^ ^ 

Br(B± ^ <i)K^) = \ ^ ' ^ (15) 

[ (9.4 ± 1.1 ± 0.7) X 10"6 (Belle) . 

The predicted branching ratios in PQCD are consistent with the experimental data. 



B. Chromomagnetic Penguin and New Physics 

The chromomagnetic penguin operator Osg in Eq. (jHl) plays an important role in the 
estimation of new physics contribution. The chromomagnetic penguin diagrams are shown 



in Fig. I3J We show only dominant diagrams in the chromomagnetic penguin amplitude 24 1 



For K meson, the moments ai and 02 were recalculated in Ref. Il 
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u, d u,d u,d 

(a) (b) (c) 

FIG. 4: chromomagnetic penguin diagrams in i? ^ (j)K decays. We show only dominant diagrams 
in the chromomagnetic penguin amphtude. 

Since there must be at least one hard gluon emitted by the spectator quark, the chromo- 
magnetic penguin amphtudes are of next-to-leading order in a<j in the PQCD formalism. 
Obviously, there are many other higher-order diagrams that must be considered simultane- 
ously. Here we limit ourselves to computing only the nonvanishing leading-order terms and 
leave the higher-order terms for future computation. As we pointed out before new physics 
contribution comes in through the chromomagnetic penguin in spite of the fact that the 
leading term is 0{a1). Of course, regular penguin amplitudes may also contain new physics 
and they are 0{as)-, and we include them. In summary, we calculate the amplitude 

A=Cl%®H{a,) + C^^,®H{as) + C^E ® H (at) . (16) 

As we will see below that these penguin amplitudes and chromomagnetic penguin amplitudes 
give comparable contributions. 

IV. NUMERICAL ANALYSIS 

In this section, we estimate the MSSM effect on the branching ratios and the CP asym- 
metries for both B*^ —>■ (f)K^ and (pK^ decays. We take single mass insertion: one 
of the LR, RL, LL, and RR insertions, which are parametrized by 6lr, Srl, Sll, and Srr, 
respectively. First, we constrain them from the branching ratio for B X^'y, which is 
an inclusive decay mode with the b ^ s transition 36]. This mode is theoretically very 
clean. The theoretical prediction with in the context of the SM agrees with experimental 
data within errors, so that this mode will give meaningful constraint on any new physics 
we might introduce. Next, we apply the constrained parameters to -B ^ (pK decays. We 
calculate the branching ratios and the direct and indirect CP asymmetries with the MSSM 
contribution. 



10 



0.02 



0.01 



-0.01 



-0.02 




E 



0.01 



0.005 



-0.005 



-0.01 





1o 


1 




2o 




[ { 


r 












a: 

■■ 0.0086 


1 1 


1 


1 



-0.04 -0.03 -0.02 -0.01 



0.005 



Re(8"LR)23 

(a) i6in)23 



-0.01 -0.005 0.005 0.01 



Re(6"RL)23 

(b) (4l)23 




^ 1 



E -1 



1 1 1 1 


1o 

2o - 






















1 1 1 1 





-1 1 



Re(8\L)23 

(C) (51^)23 



(d) Kr)23 



FIG. 5: The allowed region for {Sfj^)23, (5^/,)23> {Sll)^^^ a-iid {^rr)23 constrained from Bt{B —i- 
Xg'j). The solid lines denote the region in the la error, and the dotted lines are the 2a error. We 
take nig and rUg to be 500 GeV. In the analysis of ^ (pK, we take the following parametrization: 
(,5^^)23 = R + re'^r and {5%3 = |('^'')23|e^^ for the RL, LL, and RR insertions. 

A. Constraint from B — > Xgj 

We constrain the mass insertion parameters from Br(S — > Xg'j). The experimental result 
is Br (5 Xg'-f) = (3.3 ± 0.4) x 10~^ js^], and we take it with the 2a error to constrain the 
parameters. The results for each insertion are shown in Fig. El In this calculation, we take 
the gluino mass rrig and the squark mass rrig to be 500 GeV. The constraints for (5^/j)23 
and (5^^)23 are strong, while those for (5^2.)23 and ((5^/j)23 are very weak. This is because 
amplitudes with (5^/j)23 and (5^j^)23 are enhanced by a nig /nib factor. 
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B. B ^ <j)K Decays with the MSSM 

We calculate the MSSM effect on B° and (j)K^ decays with LR, RL, LL, 

and RR insertions. In the case of the LR insertion, we parametrize {S'[ii)23 as 

{SU),, = R + re^'^ , (17) 

where i? is a constant, and r and 9r are parameters as shown in Fig.lSfa). In the SM, (5^/2)23 
is equal to 0, that is, r = —R and 9r = 0. Since the constraints for (5^^)23 and (5^j:j)23 from 
Br{B — > Xs'y) are very weak, we take the arbitrary bounds 1(5^^)231 < 1 and | (5^^)23! < 1 
in order to use the mass insertion approximation. Both the LL and RR insertions have the 
same contribution to 5 — > (pK decays. In the case of the RL and LL{RR) insertions, we 
use the angle 9 in Fig. Efb) as a parameter: 

(5')23 = |(5')23|e^'. (18) 

In the following analysis, we scan all values on the allowed region in Fig. Efa) for the LR 
insertion. For the RL and LL{RR) insertions, we take some specific values 1(5^^)231 = 0.001 
or 0.0086, and |(5^i(^^))23| = 0.5 or 1.0. 

The Eg — Bg mixing may also be affected by the MSSM contribution '?], so that we have 
examined it in order to constraint the mass insertion parameters. The current experimental 

n 

data is AM^ > 14.4 ps~^ (at 95% C.L.) [38]. The values of AM^ is not very sensitive to 
the presence of the LR and RL insertions. In the case of the LL and RR insertions, their 
allowed regions are reduced somewhat but these insertions in the B — > (pK amplitudes is 
small so that it does not affect our analysis. 

Possible MSSM modification for B (pK branching ratios are shown in Fig.lHl As it can 
be seen from the figures, the LR insertion may give large effect on the branching ratios for 
B (f)K decays, while the contributions from RL and LL{RR) insertions are small. In the 
PQCD approach, there are large theoretical uncertainties in the calculation of the branching 
ratios. Therefore, it is difficult to obtain meaningful constraints from the branching ratios 
for B (f)K decays. In the case of the LR insertion, we suppose 16*^1 is less than 7r/2. In the 
case of the other insertions, we cannot constrain the parameters from the branching ratios. 

The results of the direct CP asymmetries are shown in Fig. [7| The experimental data of 
the direct CP asymmetry for the charged mode are A,j,K± = 0.04 ± 0.09 ± 0.01 (BaBar 3^ ) 
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FIG. 6: The branching ratios with the MSSM contribution. We take rrig and ruq to be 500 GeV. 
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and Behe, respectively. 



and 0.01 ± 0.12 ± 0.05 (Belle pal)- In the SM, the direct CP asymmetries are almost 0, 
since B (pK decays are penguin dominant processes. In the case of the LR insertion, the 
direct CP asymmetries can reach about 85%. These asymmetries are generated from the 
interference between penguin amplitudes in the SM and chromomagnetic penguin amplitudes 
in the MSSM. Since the relative strong phase between those amplitudes is large, it might be 
possible to get the large direct CP asymmetries in the PQCD approach. It is noted that the 
direct CP asymmetry in the neutral mode is almost the same as one in the charged mode. 
This result from the fact that the chromomagnetic penguin contributions as well as the SM 
contributions are almost the same in both modes. Since the current experimental data of 
the direct CP asymmetry in the charged mode is small, we expect that one in the neutral 
mode will be also small. 

Next, we consider the indirect CP asymmetry in (pKs decay. The results are shown 
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FIG. 7: Direct CP asymmetries with the MSSM contribution. We take rrig and mq to be 500 GeV. 
The dot-dashed hnes and the dotted hnes are the experimental data with the la error by BaBar 
and Belle, respectively. 

in Fig. IHl The current BaBar data is consistent with the SM prediction, however. Belle data 
is not. Our result is S(j,Ks ^ —0.28 in the LR insertion case. If we change the gluino mass 
and the squark mass, while the ratio is fixed, then S^Xg does not change. If we fix the 
squark mass and take a heavier gluino mass, then S^fiXs becomes larger. If we fix the gluino 
mass and take a heavier squark mass, then S^Ks becomes somewhat smaller, that is, new 
physics contribution becomes larger, and S^Ks remains almost constant over to be a few 
TeV. Here, we constrained (^^^2)23 from Br(i? -^s7), so that if we take a larger m^, then 
(5^^)23 increases and C^q increases slightly. In order to study the effect of the mass on 
S^iKs) "we take an extreme case where the gluino mass is 200 GeV and the squark mass is 
2 TeV. In this case, {5'Ir)2z is still (9(10~^), and the mass insertion approximation can be 
used. The result is shown in Fig. El Even in this case, the S^Ks can only reach about —40%. 
Therefore, it is difficult to explain the current Belle data by the new physics we considered 
here. 
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FIG. 8: The result of S(f)Kg with the MSSM contribution. We take rrig and rriq to be 500 GeV. The 
dot-dashed Unes and the dotted hnes are the experimental data with the la error by BaBar and 
Belle, respectively. 




FIG. 9: As an example of the extreme case, we take rrig to be 200 GeV and rriq to be 2 TeV in 
the case of the LR insertion. The dot-dashed lines and the dotted lines are the experimental data 
with the la error by BaBar and Belle, respectively. 

Finally, we comment one of the differences between our results and that of QCD factoriza- 
tion. A negative S(f,Ks in Fig-Efa) implies a negative A^Ks Fig-IZI^a). In an analysis using 
QCD factorization, a negative S^^xs implies that A^Ks{= —C^^Ks) is positive in contrast to 



our results 



. It is caused by the difference in the origin of the strong phase in PQCD and 



QCD factorization. 



V. CONCLUSION 



B — s> (j)K is one of the most important decay modes in the search for new physics. In this 
paper, we estimated MSSM contribution in the B (pK decays using the PQCD approach. 
In PQCD, strong phases are calculable, and we can predict CP asymmetries. We considered 



15 



the single mass insertions (5^/2)23, ('5rl)23) {^ll)23, and (^^/j)23, and constrain them from 
Br[B — > Xs'j). We found that the LR insertion may change the branching ratios and the CP 
asymmetries in i? ^ (pK significantly. The effect of the RL insertion is somewhat smaller 
than the LR insertion, and that of the LL and RR insertions is little. 

In the case of the LR insertion, A^i^x can reach about ±0.85 in both neutral and charged 
modes, and S^Ks ^ —0.28. The direct CP asymmetries arise from the interference between 
penguin amplitudes in the SM and chromomagnetic penguin amplitudes in the MSSM. In 
PQCD, there is a large relative strong phase between them, so that the direct CP asym- 
metries may be large depending on the new physics parameters. As in Figs. IHfa) or IHl and 
Fig. [7|^a) indicate, our result is incompatible with the current Belle data. However, the 
current Belle result is not in agreement with the current BaBar result, so that we need 
more data to arrive at a definite conclusion. Finally, it must be noted that the direct CP 
asymmetry in the neutral mode has the same tendency as one in the charged mode, because 
the chromomagnetic penguin contributions as well as the SM contributions are almost the 
same in both modes. 
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APPENDIX A: INPUT PARAMETERS 

The parameters that we used in this study are as follows js^]: \Vts\ = 0.0412, \Vtb\ = 1.0, 
\Vus\ = 0.2196, \Vub\ = 0.0036, ^3 = 80°, Mb = 5.28 GeV, Mk = 0.49 GeV, = 1.02 GeV, 
rrih = 4.8 GeV, = 174.3 GeV, Jb = 190 MeV, fx = 160 MeV, = 237 MeV, 
/J = 220 MeV, tbo = 1.54 x lO^^^ sec, tb± = 1.67 x lO^^^ g^^^ ^g)^^ ^ o.250 GeV, and 
moK = Ml/{md + m,) = 1.7 GeV. 
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